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Abstract 



- T— I 1 

^ ■ Two new sets of QCD sum rules for the nucleon axial coupling constants 



are derived using the external-field technique and generalized interpolating 
fields. An in-depth study of the predicative ability of these sum rules is 
carried out using a Monte-Carlo based uncertainty analysis. The results show 
that the standard implementation of the QCD sum rule method has only 
marginal predicative power for the nucleon axial coupling constants, as the 
relative errors are large. The errors range from approximately 50 to 100% 
compared to the nucleon mass obtained from the same method, which has 
only 10% to 25% error. The origin of the large errors is examined. Previous 
analyses of these coupling constants are based on sum rules that have poor 
OPE convergence and large continuum contributions. Preferred sum rules 
are identified and their predictions are obtained. We also investigate the 
new sum rules with an alternative treatment of the problematic transitions 
which are not exponentially suppressed in the standard treatment. The new 
treatment provides exponential suppression of their contributions relative to 
the ground state. Implications for other nucleon current matrix elements are 
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also discussed. 
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I. INTRODUCTION 



Nucleoli matrix elements of axial-vector currents at zero momentum transfer are char- 
acterized by the axial- vector coupling constants (or axial charges): isovector g^, octet g A , 
isoscalar g s A , and flavor singlet g° A . In terms of the quark spin content of the nucleon, one 
may express these coupling constants as: 

gA = Au — Ad, 

g\ = Au + Ad - 2As, 

g s A = Au + Ad, 

g A = Au + Ad + As, (1) 

where Aq is defined by {ps\q^^) 5 q\ps) = Aqu(p, s)^^y^u{p, s) and u{j>, s) is a Dirac nucleon 
spinor. Knowledge of any three of the axial charges completely determines the quark spin 
content of the nucleon. Experimentally, g& ~ 1-26 from neutron beta decay |J, g\ « 0.60 [@, 
g A ~ 0.40 from the EMC data together with strange-baryon decays |J, g A ~ 0.22 from 
SMC j| and 0.27 from E143 ||. These coupling constants reveal important information on 
the non-perturbative and non-valence structure of the nucleon. Obviously, understanding 
these quantities from QCD, the underlying theory of strong interactions, is an important 
theoretical issue. Within the framework of the QCD sum rule approach ||, these couplings 
have been previously studied in Refs. 0-0- I n ^ ne discussions to follow we will use g& in 
a generic sense to mean all of these couplings. 

Recently, a new Monte-Carlo based uncertainty analysis was introduced to quantitatively 
determine the predictive ability of QCD sum rules |15| . A comprehensive analysis of ground 



state p-meson and nucleon spectral properties was performed. Many of the findings con- 
tradicted the conventional wisdom of both practitioners and skeptics alike. In particular, 
careful consideration of opertor product expansion (OPE) convergence and ground state 
dominance revealed that the nucleon sum rule traditionally favored in the QCD sum rule 
approach is invalid. Furthermore, it was found that the nucleon interpolating field advocated 
by Ioffe is not the optimal choice. 

In this work, we will apply the Monte-Carlo method to the analysis of gA- Our goal 
is to use this new analysis tool to determine the predicative ability of the QCD sum rule 
approach for gA- This is the first application of such Monte-Carlo based analysis to a 
three-point function. Since previous studies of g& are based on the Ioffe current and the 
invalid nucleon-mass sum rule, it is important to re-investigate gA using alternate nucleon 
interpolating fields. To serve our purpose, two new sets of QCD sum rules are derived. One 
set is derived from a generalized correlator of spin-1/2 interpolating fields, and the other 
with a mixed correlator of spin-1/2 and spin-3/2 fields. We will examine all of the new sum 
rules and compare their performances. 

The organization is as follows. Sec. || sets up the basic ingredients for calculating 
using the external field method. In Sec. Ill], the QCD sum rules for the spin-1/2 correlator 
are presented and compared, wherever possible, with those obtained in previous works. In 
Sec. [TV], the QCD sum rules for the mixed correlator are given, along with a discussion of 
the phenomenological representation. Results of the Monte-Carlo based sum rule analysis 
are presented in Sec. [V|. In Sec. [VJ we give an alternative treatment of the off-diagonal 
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transitions and compare with the standard treatment. The summary and conclusions are 
given in Sec. VII . 



II. CALCULATION OF G A USING THE EXTERNAL FIELD METHOD 

The external field approach proceeds by adding an axial-vector coupling term to the 
QCD Lagrangian: 

A£ = -J2 9 q qZ^ 5 q, (2) 
q 

and considering the two-point current correlation function in the presence of a constant axial 
vector external field: 

n(p) = ijd 4 x e ip - x (Q | T{ r){x) 77(0) } I Q) z , (3) 

where \ denotes the QCD vacuum, Z is the axial vector external field and the hat notation 
denotes Z = Z a ~f a . The coupling g q keeps track of the quark flavor in the external field. 
For example, if the external current is = uj^u — d^y^d, then g u = —g d — 1, g s — 0. 
In this way, we can obtain sum rules for all the couplings in the same calculation. To first 
order in the external field, H(p) is written as: 

U(p) = U^(p) + Z a U^(p) + ---. (4) 

The axial vector coupling constant is then extracted from the sum rules for the linear re- 
sponse IlW(p). The external field couples directly to the quarks in the nucleon interpolating 
fields and polarizes the QCD vacuum. The latter can be described by the introduction of 
vacuum susceptibilities. The external field formalism is equivalent to the direct three-point 
function approach |TJJ]; the resulting sum rules are identical. 

A. Nucleon Interpolating Fields 

The most general spin-1/2 nucleon interpolating current without derivatives can be writ- 
ten as 

Vi/2 = Vi+PV2, (5) 

where j3 is a real parameter and 

Vl (x) = e abc (u aT (x)C~i b d\x)) u c {x), (6) 

rj 2 ( x ) = e abc (u aT {x)Cd\x)) l5 u c (x). (7) 

The traditional current advocated by Ioffe in QCD sum rule calculations may be recovered 
by setting /3 = — 1 and multiplying an overall factor of —2: 
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Vn = 2(77 2 - Vi) = t abc (u aT C lfl d b ) lB <fu e (8) 
However, we consider the generalized interpolator of (El) and select (3 to provide optimal sum 



rules 16 



In this work, we also consider the mixed correlator of the generalized case of Eq. (|^) : 

^,i/2 = 7m75^i/2, (9) 

with a spin-3/2 current which is also known to couple to the nucleon through its spin-1/2 
component: 

% m = e abc [ (u aT Ca pX d b )^\u c - (u aT Ca pX u b )a^d c }. (10) 

The consideration of the mixed correlator is motivated by its success in nucleon mass sum 
rules |T5| , pO| . There it is demonstrated that the sum rules from the mixed correlator have 
greater overlap with the ground state pole relative to the continuum contributions, have 
broader valid Borel regimes, and provide more stable estimates of spectral properties. It 
would be interesting to see if similar advantages can be gained here. 

The matrix elements of the currents between the nucleon and the vacuum are defined as 

(C% /2 |A0 = A 1/2 u(p) (11) 

(fik, 3/2 |iV> = A 3/2 + 7m) 7b«(p)- (12) 

where u(p) denotes the Dirac spinor of the nucleon and A describes the coupling strength of 
the currents to a nucleon state. We use the Dirac spinor normalization u{p)u{p) = 2M N . 

B. Quark Propagator in the External Field 

The calculation of the correlation function in the OPE requires the fully interacting 
quark propagator in coordinate space in the presence of the external axial-vector field. In 
the fixed-point gauge, the propagator up to order x 2 and Z is given by [0,||,|TO|,|T5 



S?(x,0;Z) = (n\T{q a (x)q b (0)}\tyz 



1 X - J) ab 

4 



X 



a a/3 + a a(3 x X nab 



2tt 2 x 

+ 32^ {9cG ^' x~ 2 2~ 

1 i xo af3 + o al3 x X nab 

+ 48 32^ [9c ap) x~ 2 2~ 

1 -(qg c a ■ Gq)a a ^ 



192 2 
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+ ^(qg c v ■ Gq) x 2 5 ab 
9g (x -Z)x-y 5 §ab 



2tt 2 

Wq 



X 1 



+ ^(qq) (xZ-x-Z) l5 S ab 
6b 



n \nab 

~^K v {qq) (Za^ + a^Z) l5 ^— 
9o 2 

+higher order terms. (13) 
The external-field-induced vacuum susceptibilities Xv an d defined by 

(qi»l,q)z = lim(z^) / Ae* 1 (Q|T 

I- H -I 

= 9 q z^Xv{qq), (14) 

and similarly, 



raft 



J29 q q( x h^5q(x) } g(o)7 M 7 5 g(o) 



\<l!lrG,nr l "<i)z = <), Z p /,'«. 6^ = ^^G 1 " 3 . (15) 



They describe the response of nonperturbative QCD vacuum to the external field. 



III. QCD SUM RULES FOR THE SPIN-1/2 CORRELATOR 

The calculation of the QCD side proceeds by contracting out the quark pairs in the 
correlation function: 

n(p) = i J d A x e ^ x (n | t{ Vl/2 (x) fj 1/a (o) } I n) z , (ie) 

where r] 1/2 is given in Eq. @, resulting in the following master formula before Fourier 
transformation: 

(Q\T{ r) 1/2 {x) fj 1/2 {0) } | n) z 

= _ e ^ c e ^' c '{ [3 2 l5 sf l5 Tr(CSf T CS b u b ' ) 

+5r'Tr(C5 cc 'JC75^' 7 5) 

+/3755r'Tr(C^ c ' T C5f 7 5) 

+{3S a u a ' l5 Tr(CS c a c ' T C l5 S bb ') 

+sf l5 cs?' T c l5 st b ' 
+Pi 5 s a / l5 csf T cs b u b ' 

+f3SfCSf T C^S bb ' 75 

+(3 2 75 SfCSf T CSt b, K }, (17) 
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where 



S?{x,0;Z) = {n\T{q a (x)q b (0)}\n} z 



(18) 



is the quark propagator given in Eq. (|13]). As discussed above, we are only concerned with 
the linear response of the correlator to the external field. When substituting the quark 
propagator into Eq. ([T7|), we keep terms only to first order in the external field. The results 
after Fourier transform have three distinct Dirac structures and can be organized by the 
invariant functions: 



z ■ n«(p 2 ) = Mp 2 ) ^75 + n 2 ( P 2 ) z ■ P p l5 + n 3 (p 2 ) iz^ 



(19) 



Three sum rules can be derived from the invariant functions. To save space, the invariant 
functions are not written out explicitly. But they can be easily inferred from the sum rules 
below. The sum rules after Borel transform are as follows. At structure Z75: 



— [([I 2 + l)(6g u + g d ) + 2/3(6^ + 5g d )} E 2 L^ 9 M 6 
64 

~ [(/3 2 + l)(9g u - g d ) + 8(3g u ] Xv a E x L~ 4//9 M 4 

\- [((3 2 + l)(6g u - g d ) + 2(3{2g u - g d )\ bE L~ 4/9 M 2 



256 
1 



+ — [{(3 2 + l){21g u + g d ) + 2/3(5^ - 3g d )] n v a E £~ 68 / 81 M 2 
4o 

[(f3 2 (18g u + 5g d ) - 2(3{2g u + 3g d ) - Ug u + g d ] a 2 L 4//9 
~2§8 [5(/3 ' + 1)9u + ^ 39u ~ 9d)] XvabL ' 4/9 



A? 



/2 



9 A (1- 



2M 



N ■ 



M 2 



+ A 



-M 2 N /M 2 



(20) 



at structure Z ■ ppjs'. 



— l({3 2 + l)(6g u + g d ) + 2(3{Qg u + hg d )\ E x L" 4 / 9 M 4 
[([I 2 + l)g d + 2P(2g u + 3g d )} XvaE L~ A ' 9 M 2 

+ 2o6 ^ + 1){69u ~ 9d) + m29u ~ 9d)] b L ~ 4/9 

[(f3 2 + l)(9g u + 5g d ) + 10(3(g u + g d )\ K v a L" 68 / 81 



144 



~ [P 2 (l8g u + g d ) + 2f3{2g u - 3g d ) + (-22g u + hg d )\ a 2 L 4 / 9 ^ 
' '? 2 + l)(2^ + 3^) + 2/5^] Xv abL- 4 / 9 ^ 



1152 



M 2 



A? 



/2 



9a . , 
M 2 



(21) 
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and at structure iZ^a^p^: 



~ [(3 2 (Q9u + 9a) ~ 2/3g d - (6g u - g d )\ a E L 2 ' 9 M 2 

+ 24 ^ 9d ~ 2/39u + 29u ~ 9 ^ Xv 0,2 

' [p 2 (6g u + 7g d ) - 2Q(3g u + 2Qg u - 7g d ] k v a 2 iT 32 / 81 — 



432 — — M 2 

' [/3 2 (2^ + 7^)-6/3^ + 4^-7^]x,^a 2 L- 14 / 27 ' 



576 1 'otj — 

-M 2 N /M 2 



A? 



/2 



g A M N 

+ A 



M 2 



(22) 



Here and in the following, a = — (27r) 2 (qq), b = (g 2 G 2 ), (qg c a ■ Gq) = —ml (qq), Xu2 = 
(27r) 2 Xi/2- As usual, the anomalous dimension corrections of the various operators are taken 
into account via the factor L = [a s (fi 2 )/a s (M 2 )] = hi(M 2 / k 2 QCD ) / ln(/i 2 /Ag CD ) , where 
fi = 500 MeV is the renormalization scale and Aqcd is the QCD scale parameter which 
will be given later. The factors E n (x) = 1 — e~ x ^Z n x n jn\ with x = w 2 /M 2 account for the 
excited state contributions, where w is an effective continuum threshold. The parameter A 
is introduced to account for all contributions from transitions between the nucleon ground 
state and the excited states; such a treatment is an approximation and may lead to errors 
in the extracted ground state property (see Sec. |IV A| below). The continuum threshold 



and transition strength are a priori unknown. So one should bear in mind that they are in 
principle different for different sum rules. We will treat them as parameters and study their 
roles in the analysis. 

The sum rules for g A can be obtained by setting g u = —g d = 1 in Eq. (p0|) to Eq. (p2|), 
while the sum rules for g A , g A , g A can be obtained by setting g u = g d = 1. Note that the 
three axial couplings g A , g A , g A share the same set of sum rules. The difference lies in the 
susceptibilities Xv and k v . 

At this point, comparisons can be made with the sum rules obtained in previous works 
using the Ioffe current by setting (3 = —1 and X\j 2 = in Eq. (p0|) to Eq. (|22|) . For the 

most part, we find that our sum rules for g A agree with those of Refs. [pj-TIDI], with only a 
few exceptions. There are differences in the anomalous dimension corrections to the terms 
involving the mixed condensate. We use —2/27 for the anomalous dimension of the mixed 
condensate, while previous works simply used 0. The dimension 7 operator Xv a b, which is 
considered in this work and in Ref. |TI|, differs at structures Zj 5 and Z ■ ppy^. Also, the 
coefficient in front of the dimension 7 operator K v a 2 at structure iZ^a^p^^ disagrees with 
that of Ref. |J. For the other axial couplings, comparisons can only be made at structure 
Z ■ pp'ys. We find that our sum rules agree with those of Refs. HJT0[|, but again report 
corrections for the XvOh term. 



IV. QCD SUM RULES FOR THE MIXED CORRELATOR 

The correlation function we consider is 
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n^(p) = i J d 4 x e^ x (Q | T{ ^, 1/2 (a;) r/„ 3/2 (0) } | ty z , (23) 

where r]^ 1/2 is given in Eq. @ and 77^,3/2 is given in Eq. fllPf). In the following, we first 
discuss the phenomenological representation of the correlation function, then calculate the 
QCD side using the OPE. 



A. Phenomenological Ansatz for the Mixed Correlator 

The linear response of the correlation function in the external field can be written as: 

Z ■ ng(p) = (-iZ a ) 1 J d 4 x e^ x (Q I T{ ^ 1/2 (x) J d A y J 5 a (y) ^ 3/2 (0) } | fl>. (24) 

After inserting two complete sets of intermediate physical states and carrying out the inte- 
grations, one has: 

ng)(p,z) = z« £ 



£<fi I r? M , 1/2 I Bps) (Bps I Jl I B'ps'}(B'ps' 1 77„ 3/2 1 0). (25) 

ss' 

The axial current coupling constant enters via the nucleon matrix element: 

(Nps I Jl I Nps) = g A u(p, s)j a j 5 u{p, s). (26) 



To determine the Dirac structure, let us look at the ground state contribution to Eq. fl25| 
after the spin sums: 



(p 2 - M 2 ) 2 

-Al/ 2 A 3 /2^ r, 2 , ,, 2 > _ A , 4^ + M 2 ,; 



/a _ M 2 f )2 L ( p + Mn ^ 7 ^ 75 Z lv + lulh^Vv - 2 7 m7 5 (Z ■ p)pj v 

c 

7^75 (Z ■ p)pp u - 2M N 7 M 7 5 (Zp - Z ■ p) j u - 8 7^75 (Zp - Z - p)p v }. (27) 



We see that there are eight distinct structures from which eight sum rules can be derived. 
The pole structure of the correlation function from Eq. ([25]) can be written as 

A1/2A3/2 9a \ - C N ^ N * x - P\- .\- , . 

(p 2 - M 2 ,) 2 + £ (p 2 - M%)(p 2 -M%.) + % (p* - M 2 *) 2 ' 1 J 

where C N ^ N * and D N ^ N * are constants. The first term is the ground state double pole, 
the second term represents the non-diagonal transitions between the nucleon and the ex- 
cited states caused by the external field, and the third term represents the excited state 
contributions. Upon Borel transform, one has 
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A1/2A3/2 9a 



M 2 



-Mff/M 2 



+ e 



-M%/M 2 



E 

.N* 



C 



N—>N* 



M 2 



N 



M 2 



N* 



1-e 



-(Ml 



-M 2 N )/M^ 



+ 



E 

N* 



M 2 N , 



(29) 



We see that the transitions (second term) give rise to a contribution that is not exponen- 
tially suppressed relative to the ground state (first term). The strength of such transitions 
at each structure is a priori unknown and is an additional source of contamination in the 
determination of g& not found in mass sum rules. The usual treatment of the transitions 
is to approximate the quantity in the square brackets by a constant phenomenological pa- 
rameter, which is to be extracted from the sum rule along with the ground state property 
of interest. Such an approximation has been adopted in the sum rules 



P to (HI). Here we 

want to stress that the transition term is in fact a complicated function of the Borel mass 
and the usual approximation alters the curvature of the phenomenological side and hence 
introduces errors in the extracted ground state property. Later in this work, we will present 
an alternative method of treating such transitions, which provides exponential suppression 
of the transitions relative to the ground state contribution. The pure excited state contri- 
butions (the last term) are exponentially suppressed relative to the ground state and can be 
modeled in the usual way by introducing a continuum model and threshold parameter. 



B. Calculation of the QCD side 

The master formula is given by 

(n\T{ % , 1/2 {x) 77„ 8/a (o) }\n) z 

= e ahc e a ' Vc '{ P^S: a ' lu a pX Tr(S b d b 'o-r> x CS c u c ' T C) 

+l P l,Sf lu a pX Tr( l5 S bb '^ x CS cc 'lc) 

-2^ 5 Sfa pX CSf T C l5 S bb '^ x 

-l,l 5 S a u a 'a pX CSf T C l5 St b 'i^ x 

-2(3 lp S a u a 'a pX CSf T CS b d b ' liy a" x 

-f3l»S a u a 'a pX CSf T CS b u b ' lu a? x }. (30) 

The calculation proceeds in the same way as in the spin- 1/2 case by substituting the quark 
propagator into Eq. (|30|), keeping terms to first order in the external field. The results after 
Fourier transform have eight distinct Dirac structures and can be organized as 

z ■ n«( P 2 ) = Mp 2 ) + Mp 2 ) wZ Pu + n 3 (p 2 ) w (z ■ vhv 

+n 4 (p 2 ) 7 M 7 5 (Z ■ p)p v + II5Q9 2 ) 7 M 7 5 ^P7v + n 6 (p 2 ) j^ZpPu 
+n 7 (p 2 ) 7 M 7 5 (Z • p)pq v + n 8 (p 2 ) 7 M 7 5 (Z • p)pp v . (31) 

After Borel transform, eight sum rules are obtained. They are as follows. At structure 

llfl^Zlv- 
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(l-(3)(g u -g d ) 



Xv aE 1 L- A ' 21 M 4 



{1 - P){9u ~ 9 -^K v aE L~ 44//81 M 2 - ^-^9u-g d ) b ^ L _ 4/27 m2 



18 



96 



. (3{7g u + 5g d ) + 10g u 2 r 20/27 , 0--P)(9*-9d) , ,-4/27 
+ g a L + XvabL 



— A1/2A3/2 



„ 2M N^ 

9A 



e -M 2 N /M 2 



(32) 



at structure / -f^ l '-f 5 Zp u : 

3(3{gu + gd) + Ug u -2g d q l8/27 m2 



18 



+ 



3/3(5fif M + 7gd) + 19^ M - 7ff d 2 r-2/9 ff(gu ~ 9d) + 7g u - 9d 2 r 8/27 



96 



+ 144 ^ m o« L — 

P(g u + llgd) - 3(8fl u + gg) ^ 2 r -8/8i 1 
+ 54 ^° 1 W 2 



A1/2A3/2 



M N K M 2 



-M 2 /M 2 



at structure 7^75 Z ■ ^7^: 

12/?(g u + g d ) + 8p u + g d L%j21 m2 _ (3(g u + 2g d ) 2 L _ 8/27 
18 ' 3 Xv 



(33) 



7(3(g u + 2g d ) 2 2 r-2/9 1 , P( 5g 9u + 55g d ) + 9(5g u - g d ) 2 r -s/8i 1 

H Z Yd Win O 7; ~r K B CI Li — 

72 A M 2 108 M 2 



A1/2A3/2 



g A M J 



N 



A 



M 2 /M 2 



e N 



M 2 

at structure 7^75 Z • pp u : 



{1 - P) [l u - 9d) Xv aE L-^M 2 



(34) 



5(1 -f3){g u -g d ) 



k*> & L 



-44/81 _ (1 ~ ^Kgu ~ 9d) fo£-4/27 



36 " 96 

4[3(2g u + g d ) + 5g u + 2<? d 

2 ^20/27 



— A1/2A3/2 



18 

9a 



M 2 



288 



Xvao m 2 



+ A 



M 2 

at structure 7^75^7^: 

3/3 (gu + gd) ~ 9u + 9d 
12 



g-M^/M 2 



(35) 



a E L 8 / 27 M 2 - + ^ ~ 7 ^ + 9d Xv a 2 L 8 ^ 27 



12 
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P(13g u + 35g d ) - 19g u + 7g d 

H Xv m n 

288 A 

/3(3l5- u + 47 g d ) - 73g u - 5g d 2 



^1/2^3/2 



216 

gAM N 

M 2 



2 a 2 L 8/27 



1 

1 

M 2 



+ A 



-Mjr/M 2 



(36) 



at structure j^sZpp^: 

(l-{3)(g u -g d ) 
24 



- 9d) + - ffrf 2 r 20/27 _L , (1- P){9u- 9d) ^ r -4/27 1 

18 a M2 + 576 Xl> 



— ^1/2^3/2 

at structure 7^75 Z • pp^ v : 

(l-[3)(g u -g d ) 
12 



4- A 
M 2+A 



e -M%/M 2 



E n L~ A ' 27 M 2 



5{l-/3)(g u -g d ) 
36 



[1- P)(9u- 9d) , j -A/27 



96 



2{l + f3) + 3g u + g d 2 20/27 1 {I - 0){g u 
9 a M 2 288 



\ v (IvLj ^ 



A1/2A3/2 



9 A 

M 2 



+ A 



-M%/M 2 



(37) 



(38) 



and at structure 7^75^ ■ ppp v : 

3(3 (g u + 9d) + llg u - 2g d r8/27 3[3(hg u + 7g d ) + l§g u -7g d 2 2/9 1 

18 96 M 2 

3P(9u + 9d) + 7g u - 9d 2 r -8/8i 1 



27 



M 4 



Ai/2^3/2 



M N M 2 



+ A 



-Mjj/M 2 



(39) 



V. MONTE CARLO SUM RULE ANALYSIS 

The reader is referred to Ref. jnjj for a complete description of the method. The ba- 
sic steps are as follows. One first generates many sets of randomly-selected, Gaussianly- 
distributed QCD parameter sets, from which an uncertainty distribution in the OPE can be 
constructed. Then a x 2 minimization is applied to the sum rule by adjusting the phenomeno- 
logical fit parameters. This is done for each QCD parameter set, resulting in distributions 
for phenomenological fit parameters, from which errors are derived. Usually, 100 such con- 
figurations are sufficient for getting stable results. We generally select 1000 which help 
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resolve more subtle correlations among the QCD parameters and the phenomenological fit 
parameters. 

The Borel window over which the two sides of a sum rule are matched is determined by the 
following two criteria: a) OPE convergence — the highest-dimension-operators contribute 
no more than 10% to the QCD side when (3 = 0, b) ground-state dominance — all excited 
state contributions (including transitions) are no more than 50% of the phenomenological 
side. Those sum rules which do not have a valid Borel window are considered unreliable and 
therefore discarded. The emphasis here is on exploring the QCD parameter space via Monte 
Carlo. The 10%-50% criteria are a reasonable choice that provide a basis for quantitative 
analysis. Reasonable alternatives to the 10%-50% criteria are automatically explored in the 
Monte-Carlo analysis, as the condensate values and the continuum threshold change in each 
sample. 



A. QCD Input Parameters 

The QCD input parameters and their uncertainty assignments are given as follows. The 
quark condensate is taken as a = 0.52 ± 0.05 GeV 3 . A number of recent studies [pl| prefer 



much larger values for the gluon condensate than early estimates of 0.47 ± 0.2 GeV 4 from 
charmonium sum rules. Hence we adopt b = 1.2 ± 0.6 GeV 4 with 50% uncertainty. The 
mixed condensate parameter is placed at m 2 , = 0.72 ± 0.08 GeV 2 . Note that the value of 
the mixed condensate itself is obtained by multiplying the randomly-selected m 2 , with the 
central value of the quark condensate. Factorization violation of the four quark operators 
is parameterized as n(qq) 2 , where we consider k — 2 ± 1 and 1 < k < 4. The QCD scale 
parameter Aqcd is restricted to the values conventionally adopted in the QCD sum rule 
approach: A QCD =0.15±0.04 GeV, and 0.10 GeV < A QCD < 0.20 GeV. Variation of A QCD 
has little effects on the results. 

The external-field-induced vacuum susceptibilities for the isovector g& have been esti- 
mated previously We consider \ v a = 0.70 ± 0.05 GeV 2 and K v a = 0.14 ± 0.14 
GeV 4 . Xv a is related to PCAC which is well known. k v is related to the matrix element 
(0\qg c G[ lu 'y u q\7(~) z which is not well determined. Note that the combinations \v a and n v a 
as a whole are selected by Monte Carlo. For the flavor singlet g A , they have been esti- 



mated in Ref. [12[]. Here we adopt them with 50% and 100% uncertainties, respectively: 
Xv a = 0.14 ± 0.07 GeV 2 , K v a = 0.01 ± 0.01 GeV 4 . 

For the octet g\ and isoscalar g s A , the susceptibilities are estimated from the following 
consideration. Nucleon matrix elements (qj^j^q) have a connected (valence) contribution, 
as well as a quark loop contribution connected only by gluons. Lattice studies [[[jj found 
that the loop contributions are almost independent of the quark flavors u, d and s. In 
the limit of equal loop contributions for u, d and s quarks, both g& and g A have the loop 
contributions cancelled out, leaving only the valence contributions, as shown in Eq. ([]]). 
Hence, we assume the susceptibilities are the same for these two couplings. On the other 
hand, g A has valence plus three flavors of loop contributions, and g s A has valence plus two 
flavors of loop contributions. From these arguments one can express the susceptibilities 
for g A in terms of those for and g° A . They are estimated as \ v = 0.33 ± 0.16 and 
k v = 0.05 ± 0.05, where we have assigned 50% and 100% uncertainties, respectively. 
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These uncertainties are assigned conservatively and in accord with the state-of-the-art in 
the literature. While some may argue that some values are better known, others may find 
that the errors are underestimated. In any event, one will learn how the uncertainties in 
the QCD parameters are mapped into uncertainties in the phenomenological fit parameters. 
Fig. [I] shows distributions for these QCD parameters drawn from a sample of 1000 sets. 



B. Search Procedure 

In principle, one can extract Qa from any of the sum rules presented earlier. In practice, 
however, some sum rules work better than others. This is because one works with a truncated 
OPE series, which may have different convergence properties at different structures. By 
selecting appropriate mixing of the components of the generalized interpolating fields, one 
can minimize the overlap with excited states, and broaden the regime in Borel space where 
both sides of the sum rule are under reasonable control. In the following, we use the optimal 
(3 = —1.2 for the spin-1/2 sum rules and (3 = for the mixed sum rules, as determined in 



Ref. [15 



The analysis of sum rules requires the corresponding nucleon mass sum rules for 
normalization. They are taken from Ref. |l5j for the spin-1/2 correlator at structure 1: 



16 16 

Note that we use the nucleon mass sum rule at the chirally-even structure 1, rather than 
the traditionally favored one at the the chirally-odd structure p, since the latter was found 
to be invalid ||15|| . The optimal fit parameters obtained from consideration of 1000 QCD 



parameter sets are 

M N = 1.17 ±0.26 GeV, \\ /2 = 0.20 ± 0.12 GeV 6 , w N = 1.53 ± 0.41 GeV, , (41) 
For the mixed correlator, we use the nucleon mass sum rules at structure 7^7^ 

\aL^E x M 4 + i±|^a6L 8 / 27 = \ 1/2 \ 3/2 M N e~W , (42) 
and at structure •y^pp u 

i a E M*- 3 -^mta - ±±j^ ^ = A 1/2 A 3/2 M N . (43) 

A combined analysis of the two sum rules from consideration of 1000 QCD parameter sets 
gives 

M N = 0.96 ± 0.08 GeV, Ai /2 A 3 /2 = 0.41 ± 0.14 GeV 6 , w N = 1.3 ± 0.2 GeV, . (44) 

Note that the parameters determined from the mixed correlator have much smaller uncer- 
tainty than from the spin-1/2 correlator. 
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FIG. 1. Distributions for the QCD parameters drawn from a sample of 1000 sets. The vacuum 
susceptibilities for are shown here. 
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When performing a Monte-Carlo analysis of a g& sum rule, we first fit the correspond- 
ing mass sum rule(s) to obtain the nucleon mass M/y, pole residue A 2 and the continuum 
threshold w^- Then, the mass and the pole residue are used in the g& sum rule where the 
three remaining parameters including the transition strength A, the continuum threshold 
w, and qa are optimized. We impose a physical constraint on w requiring that w > Mjy, 
and discard QCD parameter sets that do not satisfy this criteria. The above procedure is 
repeated for each QCD parameter set until a certain number of sets (typically 1000 sets) 
are finished. 

In doing a full search for the g^ sum rules, we encountered three scenarios regarding the 
continuum threshold w. First, w falls consistently below the nucleon mass, signaling a failure 
of the sum rule to resolve the pole from the continuum. In order to proceed in this case, the 
continuum threshold of the mass sum rule is used in the corresponding g^ sum rule. Note 
that this is the assumption made in previous works analyzing g^. Second, we are able to get 
w > Mtv from the search, but the uncertainty on w is uncharacteristically large, a sign of 
an unstable fit. The origin of the large uncertainty is that the search algorithm occasionally 
returns very large values for w. Since in the continuum model, the contributions from 
large w are exponentially damped out, the extracted results for gA and the transition are 
not seriously affected. We consider such fits marginally acceptable. One could in principle 
impose a cutoff on large w in the search algorithm to reduce the uncertainty. We choose not 
to do so, but rather use it as a performance indicator of the sum rule. Third, we get both 
w > M/v and reasonable uncertainty. This is the best scenario. 



C. Results and Discussion 

Now we are ready to examine the sum rules. We find that sum rule fl2~T|) at the chirally- 
odd structure Z ■ pfr/5, which has been chosen in previous works, fails to have a valid Borel 
window. Fig. |2| shows the highest-dimension-operator (HDO) contributions of the QCD side 
relative to the sum of terms and the continuum-plus-transition contributions relative to the 
total phenomenological side as a function of the Borel mass. The former is decreasing with 
the Borel mass while the latter is increasing. Also shown are the HDO contributions at 
(3 = 0, which should be used to determine the lower limit of the Borel window when the 
10% criteria is applied The figure reveals that around a Borel mass of 1 GeV, the 

HDO contributions are large at about 70%, indicating poor OPE convergence. Only when 
the Borel mass reaches about 1.8 GeV, do their contributions drop to below 10%. On the 
other hand, at such Borel mass the continuum contributions increase to about 80%, and 
almost completely dominate the phenomenology. These continuum model contributions are 
subtracted from the leading terms of the OPE such that the terms independent of the Borel 
mass dominate the OPE. This effect trivially explains why in Ref. [|l(]] a plateau was reached 
when the Borel mass was greater than 1.8 GeV. As a result, any ground state properties 
extracted from this sum rule are seriously contaminated by short-comings of the excited 
state model. We analyzed this sum rule using the Monte-Carlo method with 1000 sets of 
QCD parameters, and the conventional conditions: in the Borel region of 0.9 GeV to 1.2 
GeV, with the Ioffe current {(3 = —1), and the previous assumption that the continuum 
threshold in the g& sum rule is the same as that of the mass sum rule. The result is 
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0.8 1.0 1.2 1.4 1.6 1.8 
Borel Mass (GeV) 



FIG. 2. The HDO contributions of the OPE relative to the sum of terms and the contin- 
uum-plus-transition contributions relative to the total phenomenological side are displayed as a 
function of the Borel mass for sum rule @ at = -1.2 (dashed). The solid line shows the HDO 
contributions at /3 = 0. 




CD 

Oh 10 



1.2 1.4 1.6 

Borel Mass (GeV) 



1.8 



FIG. 3. The valid Borel window for sum rule (pl| ) at optimal [5 = —1.2. Both the relative 
HDO contributions limited to 10% of the OPE and the continuum-plus-transition contributions 
limited to 50% of the phenomenology are illustrated by dashed lines. Note that the lower limit is 
determined at (3 = 0. 
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gA = 11-0 ± 4.86, as compared to 1.26. This example shows the importance of maintaining 
both OPE convergence and ground state dominance in extracting hadron properties from 
QCD sum rules. 

Having demonstrated the failure of sum rule (|2l|), we now turn to sum rule (^) at the 
chirally-even structure iZ^a^p^^. This sum rule does have a valid Borel window, as shown 
in Fig. |3[ Within a wide region of about 1 GeV, the HDO contributions to the OPE are less 
than 10%, and the continuum contributions are less than 50%. Analysis of this sum rules 
yields gA = 1-87 ± 0.92, as explained below. 

In illustrating how well a sum rule works, we first cast the sum rule into the subtracted 
form lis = c gA A 2 e~ M ^l M where lis represents the OPE minus excited state contributions. 
Here c is a constant factor: 1 for @; M N for ©, (H), ©, flU, ©, 

(|36|); and 1/Mjy for fl39|). Then the logarithm of both sides is plotted against the inverse 
of M 2 . The right-hand side will appear as a straight line. The linearity of the left-hand 
side gives a good indication of OPE convergence and the quality of the continuum model. 
The two curves should match for a good sum rule. This way of matching the sum rules is 
similar to looking for a 'plateau' as a function of Borel mass in the conventional analysis, but 
has the advantage of not restricting the analysis regime in Borel space to the valid regimes 
common to both two-point and three-point correlation functions. 

Fig. |] shows the fit of ( p2|) in conjunction with ([!(]) at the optimal (3 = —1.2. Fig. |5| 



shows the fits of (p3|) and (qQ) in conjunction with (|42|) and (f43|) at the optimal (3 = 0. For 
comparison purposes, the corresponding nucleon mass sum rule is also plotted. 

Sum rule (^) does not have enough information in the OPE to completely determine 
the fit parameters, so we have assumed the equivalence of continuum thresholds in two- and 
three-point functions; w = w^. For curiosity, we also tried to fix the nucleon mass at its 
known value and found the fits and uncertainties are essentially the same. There are some 
small deviations from linearity (dashed lines) for the spin-1/2 correlator gA sum rules, but 
near perfect linearity for the mixed correlator gA sum rules. Unfortunately, the error bars 
for the Qa sum rules are much larger than those of the corresponding nucleon mass sum 
rules. 

Since all the fit parameters in the Monte-Carlo analysis are correlated, one can study 
the correlations between any two parameters by looking at their scatter plots. Fig. ^| shows 
the correlations of the QCD input parameters with gA for sum rule (|22|). The plots for 
the quark condensate a, the mixed condensate m 2 , and the vacuum susceptibility Xv a look 
fairly random, suggesting little correlation. The gluon condensate and the factorization 
violation parameter display some weak positive correlations with gA, while the vacuum 
susceptibility n v a reveals weak negative correlations. In fact, we found the same correlation 
patterns between gA, k and k v in all the sum rules. Fig. |7| shows the correlations of the 
phenomenological fit parameters with gA for sum rule (|22]). 

One should bear in mind, however, that the correlations are in general different from sum 
rule to sum rule. Fig. || shows the correlations of the QCD input parameters with for sum 
rule (PB|). Fig. |^ shows the correlations of the phenomenological fit parameters with gA for 
sum rule (|36|) . There are some subtle differences between the correlations in this sum rule 
and those in (^). First, g& is predicted to be positive here in the entire parameter space, 
unlike (|22|). Second, the distributions in Fig. || are rounded-shaped, while they are more 
concentrated toward smaller qa in Fig. |S|. Third, the correlations with the gluon condensate 
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-1 




FIG. 4. A six-parameter (-Mjv> -^i/2> W N, w, A, gX) fit of sum rule (f22^) in conjunction with (I 



The solid line is the ground state contribution: The dashed line (hidden for (40)) is the rest of the 
contributions (OPE— continuum— transitions). The error bars are only shown at the two ends for 
clarity. 




1.3 

AT 2 (GeV- 2 ) 



2.3 



FIG. 5. Same as Fig. except for the mixed correlator sum rules. From top down, the fits are 
for sum rules (p3|), (|36[), (f4"3|), and (fl^). For clarity, the mass sum rule fits are shifted downward 
by 2 units. 
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FIG. 6. Scatter plots showing correlations between g& and the QCD input parameters for sum 
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FIG. 7. Scatter plots showing correlations between and the phenomenological fit parameters 
for sum rule (p2|). 
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in (|3g) are weakly negative, opposite to those in (|22|). Fourth, the correlations among fit 
parameters are quite different in the two sum rules. It is interesting to see that in Fig. || 
small values of the pole residue are responsible for large values of qa and vice versa. 

We have examined the scatter plots for other sum rules as well. We find that the above 
differences between the mixed correlator and the spin- 1/2 correlator are qualitatively true 
for other sum rules. It is interesting to observe how different sum rules resolve the same 
observables in different ways. This reinforces the danger of linking a particular term in the 
OPE as responsible for a particular observable. 

In Table | we give a summary of the results from the consideration of 1000 sets of 
QCD parameters. Only those sum rules that have valid Borel windows are listed. The 
presence of a second row indicates that the continuum threshold was successfully searched 
as an independent parameter (the second or third scenario discussed earlier) . Note that sum 
rule (^) for g A , g A and g A has the continuum model term proportional to (3. Therefore (3 = 
—0.1 was used instead of /3 = in order to maintain a model for excited state contributions 
in place of additional poles, and enhance the magnitude of the HDO of this sum rule. 

The positive statement one can make is that for each coupling, the results are consistent 
with each other within one standard deviation. The results also agree with the experimental 
values within one and half standard deviations, although the central values appear to lie 
consistently above the experimental values. Sum rule (^) consistently predicts larger values 
for the couplings than the rest. Also the correct ordering g\ > g s A > g A is apparent. 
Unfortunately, the uncertainties in the coupling constants are large at about 50% to 100%, 
as compared to the nucleon mass obtained from the same method, which only has 10% to 
25% errors [see Eq.(f4T|) and Eq.([0])]. The large uncertainties suggest that fine adjustments 
of the QCD input parameters will allow one to make the central values of the fit parameters 
reproduce the experimental numbers. At the present stage, we feel that such refinements 
are not meaningful until the sum rules themselves are derived to a similar accuracy. 

When the continuum threshold is included as a search parameter, the only stable results 
for gA are from sum rule (0). This sum rule proves to be the best sum rule we derived. 
The overall effects of searching the continuum threshold appear to cause small increases in 
the couplings, and small decreases in the transition strengths. However, the uncertainties 
are too large to allow a definite conclusion. We have compared the quality of the fits with 
or without searching the continuum threshold and found they are essentially the same. 

The origin of the large errors lies mainly in the poorly determined normalization of the 
two-point function as indicated by the large uncertainty in A 2 . While we can consider ratios 
of two- and three-point functions to eliminate the parameter A 2 , the poorly determined 
normalization of the two-point function remains. The hope was that the two- and three- 
point functions would be sufficiently correlated that the poorly determined normalization of 
the two-point function would mirror that of the three-point function and allow an extraction 
of gA with small uncertainties. However, our analysis indicates this is not the case. 

Some interplay between gA and A 2 is illustrated in Figs. [7] and in particular §. The 
correlations exist because the valid Borel regimes for the sum rules are not identical, and 
the residue does not factor out as it would in a ratio of sum rules. Fig. || displays the 
tendency for gA to be small when A 2 is large and vice versa, as one might expect from the 
form gA A 2 c~ m nI m for the three-point function ground state contribution. We expect the 
result to be general to any QCD sum rule of nucleon current matrix elements (three-point 
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FIG. 8. Scatter plots showing correlations between g& and the QCD input parameters for sum 



rule (36) 
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FIG. 9. Scatter plots showing correlations between g& and the phenomenological fit parameters 
for sum rule (|36|). 
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TABLE I. Monte-Carlo analysis results from consideration of 1000 QCD parameter sets. A 
second row indicates that the continuum threshold could be searched as an independent parameter 
for that sum rule. Otherwise, it is assumed to be the same as that of the nucleon mass sum rule. 
The last column shows the relative contributions of the continuum (including transitions) to the 
phenomenological side in the given Borel window. 



Sum Rule 


Borel Window (GeV) 


9A 


A (GeV" 1 ) 


w (GeV) 


Continuum 


( 


22 




0.85 to 1.62 


1.87± 0.92 


0.77± 0.97 




29% to 50% 


( 


33 


) 


0.90 to 1.53 


2.52± 1.31 


0.40± 0.95 




20% to 50% 








0.90 to 1.55 


3.53± 1.44 


-0.72± 0.43 


3.49±4.73 


20% to 50% 


m 


) 


0.95 to 1.42 


1.88± 0.94 


1.20± 1.13 




31% to 50% 








0.95 to 2.34 


2.38± 1.07 


0.96± 0.92 


2.08±0.42 


20% to 50% 


Sum Rule 


Borel Window (GeV) 


9% 


A (GeV" 1 ) 


w (GeV) 


Continuum 


( 


22 


) 


0.67 to 1.16 


2.26± 2.28 


-1.49± 2.42 




22% to 50% 




) 


0.71 to 1.60 


2.55± 1.10 


0.13± 0.86 




0.6% to 50% 


PD 


1.10 to 1.23 


1.43± 0.84 


1.28± 1.06 




46% to 50% 








1.10 to 1.50 


1.57± 0.82 


1.03± 0.80 


2.97±9.25 


36% to 50% 


Sum Rule 


Borel Window (GeV) 


9% 


A (GeV" 1 ) 


w (GeV) 


Continuum 


(H 


) 


0.59 to 1.04 


1.97± 1.43 


-1.33± 0.72 




21% to 50% 




& 


0.96 to 1.95 


0.70± 0.58 


0.51± 0.63 




33% to 50% 








0.96 to 1.96 


1.06± 0.51 


0.51± 0.41 


3.72±5.40 


11% to 50% 


d 


) 


0.73 to 1.01 


0.67± 0.47 


0.94± 0.76 




34% to 50% 


Sum Rule 


Borel Window (GeV) 


9a 


A (GeV" 1 ) 


w (GeV) 


Continuum 


(H 


) 


0.64 to 0.70 


0.81± 0.99 


-1.36± 1.20 




45% to 50% 


<M 


& 


0.77 to 1.70 


0.39± 0.30 


-0.12± 0.22 




12% to 50% 








0.77 to 1.74 


0.55± 0.28 


-0.02± 0.10 


3.66± 3.15 


12% to 50% 



functions) from the linear response of the external-field. Only when the external-field is 
treated non-perturbatively, might the drawback be avoided ||14|| . 

Other sources of uncertainty are in the poorly known vacuum susceptibility k v and 
in the factorization approximation of higher order operators. From correlation studies we 
discovered that the sum rules seem to favor smaller factorization violation and larger vacuum 
susceptibility k v . Better estimates of these parameters are clearly needed. 

The presence of the transitions is an additional source of contamination not found in 
two-point functions. To see the role the transitions play, it is useful to cast the sum rule 
in the form g& + AM 2 — ■ ■ ■. By matching the two sides over some region in M 2 where a 
linear behavior is observed, one can extract g& from the intercept of the straight line with 
the y-axis. In fact this is the method explicitly used in previous analyses to extract gA, 
except that some combinations with the mass sum rule or other sum rules were used. It is 
clear that the transitions play a crucial role since they determine the slope of the line from 
which gA is extracted. Small changes in the transition strength A can lead to quite different 
values of gA- 
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VI. ALTERNATIVE TREATMENT OF THE TRANSITIONS 



As discussed in Sec. [IV A| , the transitions between the ground state and the excited 
states lead to off-diagonal contributions that are not exponentially suppressed relative to the 
ground state. We have seen in the previous section that these off-diagonal contributions play 
important roles since they determine the curvature from which gA is extracted. As can be 
seen from Eq. (|29|) , the contribution from the transitions is in general a complicated function 
of the Borel mass. The usual treatment of crudely modeling the transitions by a constant 
parameter (multiplied by e~ M ^/ M ) could have significant impact on the determination of 
qa- Recently, a new formalism for treating the transitions has been pointed out in Ref. [E2 



In this new formalism, the transition contribution is exponentially suppressed relative to 
the ground state contribution and hence can be included in the continuum model. Here we 



briefly illustrate the new formalism. The reader is referred to Ref. [f22[ for further details 
and discussions. 

Let us consider the phenomenological representation of Eq. (|28|). If one first multiplies 
the expression with the factor (p 2 — Mjy), then performs the Borel transform, one sees 
that the transition contributions are exponentially suppressed relative the ground state 
contributions. As a result, one can use the conventional pole plus continuum model on 
the phenomenological side. The physics behind such a manipulation is a rearrangement 
of information: the transition strength on the phenomenological side is not lost, but gets 
absorbed into the new continuum model built from an altered OPE. So it is important in 
this new formalism to treat the continuum threshold as an independent phenomenological 
parameter to be extracted from the sum rule along with the ground state property of interest. 
This point has been emphasized in Ref. p2" |, where an example was given in the scalar 



channel. (This point, however, was completely ignored in Refs. p3|-p5|, where a similar 
formalism was adopted). 

The sum rules in this new formalism can be obtained from the standard ones by the 
following substitutions. On the OPE side, M 6 -> 3M 8 - M^M 6 , M 4 2M 6 - M 2 N M A } 
M 2 -> M 4 — MfjM 2 , M° -> -M 2 N , 1/M 2 -> -(1 + M^/M 2 ), 1/M 4 -> -(1/M 2 + M^/2M 4 ). 
The E factors associated with powers of the Borel mass are suppressed for clarity. On 
the phenomenological side, for the spin-1/2 correlator, one replaces the right hand side of 
Eq. (^Up to Eq. (^) by c X 2 1/2 g A e~ M H M \ where c is 2M 2 N , -1, -M N , respectively. Similarly 

the phenomenological side of the mixed correlator becomes c A1/2A3/2 gA e~ M x/ M2 , where c is 
2M 2 N} 2M N , -M N , -1, -M N -1,-1 -1/Mjv, for Eq. (H) to Eq. (g|), respectively. We will 
denote these modified sum rules by appending a symbol (PC) indicating pole+continuum 
phenomenology to their standard counterparts. 

In Table [TJ we summarize all the results from a Monte-Carlo analysis of 1000 QCD 
parameter sets using the new formalism. For comparison purposes, we also list the results 
from not searching the continuum threshold. 

The results for the coupling constants are qualitatively the same as those in the standard 
approach. This suggests that the main source of error in is not in the treatment of transi- 
tions. However, the new formalism does show the potential for improvement. For example, 
the Borel windows become generally wider, the continuum contributions become smaller, 
and the numerical fits are generally more stable. Future studies should take advantage of 
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TABLE II. Same as Table |, except they are obtained from the new pole-plus-continuum sum 
rules, as indicated by (PC). 



Sum Rule 


Borel Window (GeV) 


9A 


w (GeV) 


Continuum 


©(PC) 


0.82 to 1.83 


1.44± 0.59 




2.3% to 50% 


©(PC) 


0.84 to 1.40 


2.80± 1.50 




9% to 50% 




0.84 to 1.65 


3.60± 1.76 


2.22± 0.62 


0.5% to 50% 


©(PC) 


0.95 to 1.60 


2.02± 1.02 




9.3% to 50% 




0.95 to 1.81 


2.44± 1.12 


2.15± 0.48 


10% to 50% 


Sum Rule 


Borel Window (GeV) 


9 8 A 


w (GeV) 


Continuum 


©(PC) 


0.65 to 1.36 


1.40± 1.03 




2% to 50% 




0.65 to 1.65 


1.75± 1.22 


2.48± 2.08 


0% to 50% 


©(PC) 


0.67 to 1.48 


2.44± 1.30 




2% to 50% 




0.67 to 1.77 


3.19± 1.54 


2.39± 0.57 


0% to 50% 


©(PC) 


1.13 to 1.96 


1.43± 0.75 




9% to 50% 




1.13 to 2.17 


1.65± 0.84 


2.07± 0.78 


9% to 50% 


Sum Rule 


Borel Window (GeV) 


9l 


w (GeV) 


Continuum 


©(PC) 


0.52 to 1.09 


0.62± 1.48 




1% to 50% 




0.52 to 1.33 


0.54± 0.99 


2.16± 2.26 


0% to 50% 


©(PC) 


0.57 to 1.24 


1.20± 0.94 




2% to 50% 




0.57 to 1.50 


1.56± 1.08 


2.39± 2.15 


0% to 50% 


©(PC) 


0.98 to 1.68 


0.79± 0.63 




9% to 50% 




98 to 1 90 

W • xj K-) \J\J -L • > J \ > 


97± 62 


2 30± 97 


1% to 50% 


Sum Rule 


Borel Window (GeV) 


9a 


w (GeV) 


Continuum 


©(PC) 


0.45 to 1.03 


-0.36db 0.90 




0.3% to 50% 


©(PC) 


0.60 to 0.91 


0.36± 0.59 




9% to 50% 




0.60 to 1.11 


0.56± 0.59 


2.22± 1.79 


0.5% to 50% 


©(PC) 


0.76 to 1.44 


0.42± 0.32 




5% to 50% 




0.76 to 1.68 


0.51± 0.36 


2.16± 1.27 


5% to 50% 



the new formalism to avoid the errors associated with the approximation in modeling the 
unwanted transitions. 

Although the extracted couplings appear to agree within uncertainties with or without 
searching the continuum thresholds, there are qualitative differences between the two. To 
demonstrate this, we show in Fig. |l0]the fits for sum rule (|22])(PC), (|33f)(PC) and (|36|) (PC) , 
without searching the continuum threshold. In Fig. [11] we show the two successful fits 
obtained when searching the continuum threshold. We see that the two sides of the sum rules 
do not match very well when the continuum threshold is not searched. The match improves 
significantly for the mixed correlator sum rules when the continuum threshold is searched. 
A large uncertainty was found for the searched continuum threshold in sum rule fl2"2])(PC). 
This suggests that it is really unable to separate the pole from the continuum. 

The above results confirm that the continuum threshold in the new formalism is relied 
upon to compensate for the transition strength that was explicitly taken into account in 
the standard approach. Fixing it in the new approach would introduce an artificial bias 
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FIG. 10. From top down, Monte-Carlo based fits for gA from sum rules (|33[)(PC), (|36|)(PC), 
and (|2^)(PC). The continuum thresholds are assumed the same as that of the corresponding mass 
sum rules. The solid line is the ground state contribution, and the dashed line is the OPE mi- 
nus continuum. For clarity, the fits for (|36|)(PC) and (|2^)(PC) are shifted downward by 1 unit, 
respectively. 




FIG. 11. From top down, fits of sum rules (|33|)(PC) and (|36|)(PC) in which the continuum 
threshold was searched successfully. 
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to the extracted ground state spectral properties. This is in contrast with the situation in 
the standard approach. There the OPE was not altered and the transitions were modeled 
explicitly. As a result, the quality of the fits with or without searching the continuum 
thresholds was essentially the same. 



VII. SUMMARY AND CONCLUSIONS 

We have derived eleven new QCD sum rules for the nucleon axial vector coupling con- 
stants using the external-field method and generalized interpolating fields. Three are from 
the spin-1/2 correlator, which can be reduced and compared to those in the literature. Eight 
are from the mixed correlator, which are new. Using the Monte-Carlo analysis, we are able 
to determine the predicative ability of the new sum rules for for the first time. 

The main advantage of the Monte-Carlo analysis is that it takes into account all un- 
certainties in the QCD input parameters simultaneously, allowing a quantitative study of 
how the uncertainties in the QCD input parameters propagate to the phenomenological fit 
parameters. The method carefully monitors the OPE convergence and the ground-state 
dominance, the two key criteria in order for the QCD sum rule method to work. Together 
they determine the Borel window over which the two sides of the sum rules are matched. 
Those sum rules which do not have a valid Borel window are considered unreliable and 
therefore discarded. Conventional QCD sum rule analysis is limited to only a small portion 
of the QCD parameter space, and often the uncertainties assigned to the fit parameters are 
not based on rigorous error analysis but rather have a certain degree of arbitrariness. 

Our most important findings in this work are: a) The nucleon axial vector coupling con- 
stants calculated from standard QCD sum rule method have large uncertainties associated 
with them, approximately 50% to 100%, as compared to the nucleon mass obtained from 
the same method which has only 10% to 25% error. Within uncertainties, the numbers 
obtained from different sum rules are consistent with each other and with the experimental 
values at 1.5er. The correct ordering in magnitude of the coupling constants, g\ > g s A > g% 
is also predicted by the sum rules, b) Sum rule (0), upon which the previous analyses 
of qa are based, has poor OPE convergence and poor ground-state dominance properties. 
The results extracted from this sum rule are unreliable. Both of these findings contradict 
the conventional wisdom. Traditionally, 10% to 20% errors are often claimed for g& from 
QCD sum rule calculations without rigorous error analysis. The selection of sum rule ( |2T|) 
was based upon similarities with the mass sum rule M at the structure p, which was shown 



to have poor convergence properties fL5[ , or dimension arguments [Tfl. Our Monte-Carlo 
analysis has shown that such criteria may not be reliable in selecting sum rules. 

The origin of the large errors in g& is three- fold. First, it is in the poorly determined 
nucleon mass sum rules, which are used to normalize the couplings extracted from the form 
gA A 2 6~ m n/ m . Second, two new parameters are needed to describe the response of the QCD 
vacuum to the external field. These new parameters introduce additional uncertainties, es- 
pecially the vacuum susceptibility k v which is very poorly known. Third, the transitions 
between the ground state and the excited states caused by the probing axial current is an- 
other source of uncertainty. Their contributions are not exponentially suppressed relative to 
the ground state and must be included in the spectral representation. Little is known about 
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these transitions. The standard approach is to introduce a new phenomenological parame- 
ter to account for all the contributions from the transitions. Usually the new parameter is 
assumed to be constant, but in fact it depends in some complicated way on the Borel mass. 
The approximation may have a sizable impact on the results since it plays a crucial role in 
the extraction of g&- In any event, the presence of an additional unknown parameter leads 
to larger uncertainties in g^. 

To further investigate the role of the transitions, we have studied an alternative treat- 
ment which can lead to exponential suppression of their contributions relative to the ground 
state. As a result, one can apply the traditional pole-plus-continuum model in the phe- 
nomenological representation. The contributions of the transitions are not lost, but simply 
get absorbed in the new continuum model. So it is important in this formalism to search the 
continuum threshold as an independent parameter. The results for the coupling constants 
are similar to those in the standard approach. This implies that the main source of error is 
not in the treatment of transitions, but in the pole residues and the vacuum susceptibilities. 
However, there are advantages to be gained with the new formalism. The Borel windows 
become generally wider, and the numerical fits are generally more stable. Future studies 
should take advantage of the method to avoid errors associated with the approximation of 
modeling of the unwanted transitions. 

As for the question of valid sum rules for g^, our Monte- Carlo analysis reveals that only 
three out of the eleven sum rules are able to resolve the ground state properties from the 
continuum. They are sum rule (p2|) from the spin-1/2 correlator, (|33|) and (|36|) from the 
mixed correlator. Even the three valid sum rules perform differently. Sum rule fl36|) has 
the best performance in terms of allowing a full search and the stability of the results. It is 
followed by sum rule (|33|) , then by (p2|). We find that in general the mixed correlator sum 
rules perform better than the spin-1/2 correlator sum rules, similar to the situation for the 
nucleon mass sum rules. 

Study of correlations among the input and fit parameters reveals how a particular sum 
rule resolves the ground state from the continuum. We find that gA does not have strong 
correlations with the quark condensate, the mixed condensate, the gluon condensate, nor 
the vacuum susceptibility Xv 9a has some weak positive correlations with the factorization 
violation parameter k, and some weak negative correlations with the vacuum susceptibility 
k v . These correlations can give some hints on what values for k and k v are preferred by the 
self-consistency requirement in the sum rules. We should stress that correlation patterns are 
different for different sum rules. The above observations are some general trends displayed 
by all the valid sum rules, despite some subtle differences. 

As to the problem of how to reduce the large uncertainties in g^ from QCD sum rule 
calculations, one may attack it from several directions as suggested from the discussion of 
their origins: a) improved accuracy of the QCD sum rules, b) better estimates of the vacuum 
susceptibilities, especially k v , and c) better treatment of the transitions. Of these, point a) 
is the most important. One way to improve the accuracy of the sum rules is to reduce the 
uncertainties in the QCD input parameters. However, there is a limit on it. As shown in 



Ref . |15| , the best one can do with current implementation of the QCD sum rule approach is 
25% for the two-point function normalization. Reduction of this uncertainty to the 10% level 
requires 5% uncertainties on the QCD input parameters. But this level of accuracy is beyond 
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the current state of the art as indicated by the unacceptably large xV^DF - We believe the 
future lies in improving the sum rules themselves. Examples may include incorporating a s 
corrections, higher order power corrections, and devising alternative interpolating fields that 
result in better sum rules, as done here. 

Finally, we want to point out that although this work is focused on qa, the conclusions 
are expected to apply to other nucleon current matrix elements. The presence of large 
uncertainties in nucleon matrix elements indicate that only small adjustments of the QCD 
input parameters are required to reproduce the measured values. Likewise, small fluctuations 
in the input parameters can lead to large fluctuations in the matrix elements, such that 
most QCD sum rule calculations of matrix elements in the literature are likely to have 50% 
uncertainties at best, given the current implementation of the QCD sum rule formalism. 
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